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We demonstrate that the flow of a longitudinal spin current with different spin polarization 
will induce different patterns of charge accumulation in a two-terminal strip, or electric current 
distribution in a four-terminal Hall-bar structure, of two-dimensional electron gas with Rashba spin- 
orbit coupling (RSOC). For an in-plane polarized spin current, charges will accumulate either by the 
two lateral edges or around the center of the strip structure while, for an out-of-plain polarized spin 
current, charge densities will show opposite signs by the two lateral edges leading to a Hall voltage. 
Our calculation offers a new route to experimentally detect or differentiate pure spin currents with 
various spin polarization. 

PACS numbers: 85.75.-d, 72.20.My, 71.10.Ca 



Semiconductor spintronics has achieved remarkable 
success in the past decade and is still progressing rapidly. 
And spin-orbit science and engineering, which allow for 
electrical manipulation of spin polarization and spin cur- 
rents in nonmagnetic semiconductors, is one of the key 
steps to implement spintronic devices Utilizing spin- 
orbit coupling, various schemes were proposed to gener- 
ate pure spin currents,- while to detect pure spin currents 
remains a challenge from either experimental or theoret- 
ical aspect. The detection of spin currents often involves 
spin accumulation or electrical effects, despite quantum 
interferences by optical means also reported.'^ For exam- 
ple, spin accumulation induced by spin current near the 
boundaries has been detected in both n- and p-doped 
semiconductor systems experimentally,'* and Hall voltage 
resulted from the reciprocal extrinsic spin Hall effect was 
observed in diffusive metallic conductors.-'- Recently it 
was reported that an optically injected spin current flow- 
ing through a Hall-bar semiconductor can generate either 
in- ward or out-ward electric currents while the Hall volt- 
age remains zero J- This observation renders a manifesta- 
tion of the tensor-like nature of the spin current, which 
means that both the spin polarization and the direction 
of motion are decisive factors in producing physically ob- 
servable effects. 

The electrical detection of spin currents has reduced 
complexities in practice and thus is potentially more 
applicable i^i^ Whereas a fundamental problem of theo- 
retically studying the electrical effects resulted from spin 
currents, is the difficulties to incorporate the concept 
of spin current into a theoretical formalism, for reasons 
like the ambiguous definition of spin currents under cer- 
tain circumstances;^ In this Letter we investigate such 
effects in a mesoscopic system of spin-orbit-coupled two- 
dimensional gas (2DEG) with ideal leads which connect 
several special electron reservoirs. Generation of the spin 
current is simulated phenomenologically by introducing 
spin-dependent chemical potentials within each electron 



reservoir. These chemical potentials are tuned separately 
for each spin component to produce independent poten- 
tial gradients, so that electrons with opposite spins are 
driven to move in opposite directions through the spin- 
conserved leads and the spin-orbit-coupled central region. 
It must be emphasized that the aim of this phenomeno- 
logical simulation is to capture only the general and es- 
sential features of the spin current flow, without bother- 
ing about the details of the method in its generation and 
injection, and this is justified as long as we keep our fo- 
cus on the electrical effects that are led by or intimately 
related to the circulation of the spin current. 

We demonstrate spin-current-induced electrical effects 
by showing charge accumulation induced in a strip struc- 
ture and the electric current distribution in a correspond- 
ing Hall-bar structure. The Landauer-Biittiker-Keldysh 
formalism is used in our calculations, which is a quantum- 
in-nature approach and has extensive applicationi^ The 
total Hamiltonian of the central region is He — Hg+Hso 
, where Hq = ■|;^+Vb is the kinetic energy plus the hard- 
wall confining potential Vq, m* is the effective electron 
mass, and Hso = —a{k x cr) ■ z, with a the strength 
of Rashba spin-orbit coupling (RSOC), k the wave vec- 
tor, (T the vector composed of Pauli's matrices, and z 
the unit vector perpendicular to the plane of the 2DEG. 
After discretizing He with the tight-binding approxima- 
tion and transforming it into the spin bases which are 
the eigenstates (denoted by Xm O'" Xi' where /i, v =t or 
I) of f • cr, where f is the orientation of spin polarization 
under consideration, we have 

<ij> Ai," 

where and c are the creation and annihilation opera- 
tors of electrons at sites i(j) with spin ^{v) and < ■ • • > 
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in the central region, and km is the wave vector along the 
semi-infinite lead. Here we assume that there is no spin- 
orbit coupling in the leads. This not only guarantees that 
the spin currents under our investigation is well-defined 
from experimental aspect, but also is justified because 
it turns out that the interface mismatch between the 
leads and the central part actually contributes little to 
the patterns we observed. The retarded Green's function 
C'iE) ^{E-Hc~ and the lesser Green's func- 

tion is given by the Keldysh equation — G"'E<G°, 
with E<{E) = -Ep.^fiE - e^^)mjE) - ^JE)), 
where is the spin-dependent chemical potential for 
electrons of spin ^ in the lead p, and f{E — e!^) is the 
Fermi-Dirac distribution function. Expressed in terms 
of the lesser Green's function, when a steady state is 
achieved, the non-equilibrium charge density at site i is 

/oo 
dEG<{i,^l■^,^i;E) (3) 

and the electric bond current from site i to site j is 
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FIG. 1: Spin-current-induced charge distributions in a 
200nmxl00nm strip, where the spin current is polarized 
along y and the distributions are displayed after the spin- 
independent background distribution (when a = 0) has been 
subtracted from, (a) the image of charge distribution when 
a = —6.1 X 10~^^eVm. (b) the image of charge distribution 
when a = 6.1 x 10~^^eVm. (c) and (d) the averaged (in terms 
of the longitudinal dimension) charge distributions along the 
narrow side of the sample, comparing different cases with a 
having the same sign (negative for plot (c) and positive for 
plot (d)) but increasing magnitudes. 



means the pairs of nearest neighboring (n.n.) sites, and 



''ij 



i = j ± 4 



where ^^(j,) is the unit vector displacement between two 
n.n. sites in x{y) direction, to = /2m*a? is the over- 
lap integral of two n.n. sites with a the average spacing 
between two n.n. sites, t^Q = a/2a, / is the identity 
matrix, and u = (xf , Xl) is the unitary matrix which ro- 
tates az to f -cr. To count in the effect of the semi-infinite 
ideal leads, self-energy terms are introduced, Y7 = 
with the specific term due to lead p in the n spin diagonal 
block is 



where (pm{Pi) is the m*^ eigenfunction in the transverse 
dimension at site pi in lead p which is adjacent to site i 



Our present study focuses on electrical effects resulted 
from spin current with spin polarization orientated in 
plane (referred to as y henceforth) which is perpendic- 
ular to the direction of motion (referred to as x hence- 
forth) of electron spin. When a spin current with y- 
spin-polarization is flowing through a strip, depending 
upon the sign of the RSOC a, which stands for whether 
the velocity, the spin polarization and the gradient of 
potential for the RSOC satisfy the left-hand or right- 
hand chirality, the electric charges will accumulate either 
near the lateral edges or around the middle of stripe, 
as shown in Fig[TJ In this figure, two equal-magnitude 
-f y/— y-spin-polarized current are driven in the direc- 
tion +x/~-x by the spin-dependent chemical potentials 
e^.^ ~ ± sign(/i)0.1io through a 40 x 20 lattice with the 
Fermi Energy Ef = O.lio, which is measured from the 
bottom of the conduction band and is small enough to 
ensure the parabolic energy-momentum dispersion in the 
tight-binding approximation, and a — —6.1 x 10~^^eVm 
or t|o = -0.02io in FigH^a), and a = +6.1 x lO'^^eVm 
or tgQ — +0.02to in FiglTJb), respectively. It is easy to 
identify that the charges of carriers accumulate by the 
two lateral edges in the former case and at the middle 
part in the latter one. This is further manifested by a 
comparison of Fig[TJc) and Figlljd), where the averaged 
charge densities with ±a are plotted. For a negative a, 
increased magnitude will lead to increased accumulation 
by the edges, while for a positive a, the same trend hap- 
pens at the middle. It should be noted that since we 
do not consider any dissipation mechanisms inside our 
samples, these accumulations should be interpreted as 
an effect purely due to quantum coherent transport. 
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FIG. 2: Charge current distributions after two additional 
leads are connected to the lateral edges of the strip in Fig. 1. 
The direction and magnitude of local current density are in- 
dicated by the orientation and length of arrows, respectively. 
The transverse electric currents are both flowing outward in 
plot (a), where a — —6.1 x 10~^^eVm and the ratio of the 
total induced electric current to the circulating spin current 
Ic/Is = 2.57 X 10"^ 
where a = 6.1 x 10~ 



-j-p^, or both flowing inward in plot (b), 
'^^eVm and h/h = 1.14 x 10"^^. 
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To show physical consequences of charge accumulation 
in a strip structure, we calculated electric current dis- 
tributions in a Hall-bar structure where two additional 
leads are symmetrically attached to the two lateral sides 
of the original strip as illustrated in FigO The current 
distribution for either sign of a in this structure has re- 
markable consistence with the charge accumulation in 
the corresponding strip system in terms of that the ad- 
ditional leads act just as the pathways for the accumu- 
lated charges to flow through. Specifically, when a < 0, 
charges accumulated by the lateral edges will flow out- 
ward through the two transverse leads and at the same 
time the net charges will be drawn though the two longi- 
tudinal leads, while when a > 0, the charges tend to be 
drawn inward through two transverse leads and flow out 
through two longitudinal ones, which accounts for the 
accumulation of charges around the center of the strip- 
shaped sample. Quantitatively, the induced electric cur- 
rent in each transverse lead sums up to be typically two 
orders less than the magnitude of the total circulating 
spin current. These observations are also fully consistent 
with the electric current patterns in referenceiS where 
the linear response approximation are adopted to pro- 
duce the results. 

So far we have mainly discussed the charge density and 
the electric current distribution induced by a spin current 
with in-plane spin polarization (along y). Yet a tensor 
as a spin current is in essence, it is also worthwhile to 



FIG. 3: Spin-current-induced charge distributions in a 
200nmx lOOnm strip, where the spin current is polarized along 
X in (a) and along z in (b). a = 6.1 x 10~^^eVm and the distri- 
butions are displayed after the spin-independent background 
distribution has been subtracted from. 



investigate the spin current with different configurations 
of spin polarization and velocity. Without losing the gen- 
erality, we concentrate on three cases of spin along x, y 
and z, respectively. From the study of these special cases, 
the generic result for a spin current with arbitrary spin 
polarization can be derived, and also maximum symme- 
tries can be observed therein. Since one of these cases 
with y-spin-polarization has been presented already, we 
show the charge distributions induced in other two cases 
in Fig [3) Compared with the induced charge distribu- 
tion shown in FiglTJb), which is calculated with all the 
same parameters except for the spin polarization of the 
spin current, the induced charge distributions in Figl3]Ja) 
and (b) show clear differences in terms of the spacial 
symmetries they possess, that is the C2v symmetry in 
the case when spin polarized along y, the C2 symme- 
try when spin polarized along i, and the Cs symmetry 
when spin polarized along This is because the un- 
derlying system is invariant under each operation of the 
C2V space group combined with an appropriate unitary 
transformation in the spin space, while in order that the 
circulating spin current is also invariant under that com- 
bination of transformations, the valid space-group oper- 
ations will be limited into a subgroup of C2V (C2u it- 
self in the y-spin-polarization case, and C2 or Cs in the 
X- or i-spin-polarization case respectively), accordingly 
the charge distributions induced by different spin cur- 
rents will exhibit different symmetry properties. Also it 
is easy to infer that the electric current distributions or 
any other electric effect induced by a corresponding spin 
current will display the same spacial symmetry as long 
as the structural symmetry of the system is preserved, 
which is actually a general yet rigorous limitation and 
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has been justified in all our calculations that are shown 
or not shown here. 

Moreover, the relevance of the spin polarization of spin 
current to the induced charge accumulation also lies in 
the differences between the average distributions along 
the transverse of the strip resulted from spin currents 
which are polarized in different directions. In particular, 
Figdjc) and (d) show that the charges of carriers tend to 
accumulate, depending on the sign of a, either by both of 
the two lateral edges or around the center of the sample 
when the spin current is polarized along y. When the 
spin current is polarized along z, however, it can be ob- 
served from FiglSlJb) that the charge accumulation on av- 
erage will have opposite signs near the two lateral edges, 
which is closely related to as well as consistent with the 
reciprocal version of the spin Hall effect, and was re- 
ported to be observed in platinum wires. ^ In another case 
that the spin current with i-spin-polarization is consid- 
ered, the pattern of the charge accumulation (shown in 
FigEIa)) is more complex in the sense that there are re- 
versed accumulating trends at the two sides of the trans- 
verse plane passing through the symmetric center, which 
is determined by the C2 symmetry of this system men- 
tioned above. While after average has been taken over 
the longitudinal dimension, it can be reasonably expected 
that the opposite accumulations of two halves will tend 
to cancel each other and produce a weakened net effect 
with symmetrical distribution about the transversal cen- 
ter. In contrast to that in the y-spin-polarized case, in ei- 
ther case with the x- or £-spin-polarization, the averaged 
distribution along the traverse does not change under the 
reversal of the sign of a. Besides, we notice that the elec- 
tric effects induced by different spin currents are actually 
contributed by different ranks of power with respect to a, 
which is also manifested when the sign of a is reversed. 
Specifically, when the spin current is polarized along x or 



y, it is mainly the linear a that is responsible for the in- 
duced electric effects, while in the z-spin-polarized case, 
it is the second rank, i.e. , playing the role. 

We conclude this letter by evaluating the accessibility 
of an experimental observation to the electrical patterns 
investigated here. From the data plotted in Figlljc), 
we estimate roughly the electrostatic potential differ- 
ence Ay to be of order O.lmV between either of the 
peaks and the valley of the charge distribution as long 
as a > 3 X 10~^^eVm. And in a realistic semiconduc- 
tor quantum well with the Fermi energy typically being 
tens of meV, the electronic potential difference will in- 
crease to at least several meV or tens of Kelvins subject 
to increased RSOC strength. This implies that the tem- 
perature requirement for observing these patterns can 
be well satisfied within current experimental conditions. 
Regarding these features as well as the sample size, we 
propose the use of Kelvin probe force microscopy^'^ in ob- 
serving the charge accumulation predicted here. On the 
other hand, we point out that the present work may also 
account for a recent experimental observation of the spin- 
current-induced electric currentsji which possess the key 
features exhibited in Fig[T]and Figj^l that is, there is no 
Hall voltage while the electric currents circulate through 
x-channels to y-channels. And quantitatively the experi- 
ment is consistent with the calculated ratio of the induced 
charge current to the spin current captioned in Fig[TJ In 
short, the spin current with in-plane spin polarization 
may produce measurable charge accumulations or elec- 
tric currents with a novel behavior, and the underlying 
effects may open a promising way to the electrical detec- 
tion of spin currents. 
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